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Abstract 

The majority of recently emerging infectious diseases in humans is due to cross-species pathogen transmissions from 
animals. To establish a productive infection in new host species, viruses must overcome barriers to replication mediated by 
diverse and rapidly evolving host restriction factors such as protein kinase R (PKR). Many viral antagonists of these restriction 
factors are species specific. For example, the rhesus cytomegalovirus PKR antagonist, RhTRSI, inhibits PKR in some African 
green monkey (AGM) cells, but does not inhibit human or rhesus macaque PKR. To model the evolutionary changes 
necessary for cross-species transmission, we generated a recombinant vaccinia virus that expresses RhTRSI in a strain that 
lacks PKR inhibitors E3L and K3L (VVAEAK+RhTRSl). Serially passaging VVAEAK+RhTRSI in minimally-permissive AGM cells 
increased viral replication 10- to 100-fold. Notably, adaptation in these AGM cells also improved virus replication 1000- to 
10,000-fold in human and rhesus cells. Genetic analyses including deep sequencing revealed amplification of the rhtrsl 
locus in the adapted viruses. Supplying additional rhtrsl in trans confirmed that amplification alone was sufficient to 
improve VVAEAK+RhTRSI replication. Viruses with amplified rhtrsl completely blocked AGM PKR, but only partially blocked 
human PKR, consistent with the replication properties of these viruses in AGM and human cells. Finally, in contrast to AGM- 
adapted viruses, which could be serially propagated in human cells, WAEAK+RhTRSI yielded no progeny virus after only 
three passages in human cells. Thus, rhtrsl amplification in a minimally permissive intermediate host was a necessary step, 
enabling expansion of the virus range to previously nonpermissive hosts. These data support the hypothesis that 
amplification of a weak viral antagonist may be a general evolutionary mechanism to permit replication in otherwise 
resistant host species, providing a molecular foothold that could enable further adaptations necessary for efficient 
replication in the new host. 
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Introduction 

There are at least 868 described zoonotic microbial pathogens, 
33% of which are capable of human to human transmission [1]. 
Recent viral zoonoses have led to some of the most devastating 
and medically relevant outbreaks in modern history, including 
SARS coronavirus, pandemic influenza, and HIV/AIDS, high- 
lighting the urgent need to understand how viruses adapt to infect 
new species. At a population level, factors influencing the 
transmission of zoonotic pathogens to humans include increasing 
population density, greater contact with wildlife, increased travel, 
and poor public health infrastructure [2,3] . However, these factors 
only allow the microbe increased access to new hosts; they do not 
directly enable it to adapt to and replicate in the new species. 
Intermediate hosts, animals that are not the natural host of a virus 
but are still permissive or semi-permissive for viral replication, play 
a critical role in cross-species transmission. These hosts can 
facilitate increased contact between a virus and a new host, and 
drive adaptive changes that may improve virus replication 
(Reviewed in [4]). For example, spill-over of Nipah virus from 



fruit bats into pigs, the intermediate host, increased human 
exposure to the virus and resulted in eventual human outbreaks in 
Malaysia [5,6]. In another example, lentiviral adaptation through 
intermediate chimpanzee hosts led to both increased contact with 
humans, and adaptive genetic changes permitting the virus to 
inhibit the human versions of several host restriction factors 
(Reviewed in [7]). 

At a molecular level, the initial success of a virus after entry into 
a new host cell depends on its ability to overcome cellular host 
restriction factors. A subset of these proteins inhibits specific virus 
families, such as the restriction of retroviruses mediated by 
TRIM5a [8] . However, other restriction factors, including protein 
kinase R (PKR), block the replication of multiple virus families. 
PKR is activated by binding to double-stranded RNA (dsRNA), a 
common byproduct of both RNA and DNA virus replication, 
followed by dimerization and autophosphorylation. Activated 
PKR then phosphorylates the Gt-subunit of eukaryotic initiation 
factor 2 (eIF2a), ultimately arresting translation initiation [9]. In 
response to the broad and potent barrier to viral infection imposed 
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Author Summary 

The spread of microbes from animals to humans has been 
responsible for most recently emerging human infectious 
diseases, including AIDS, bird flu, and SARS. Therefore, 
understanding the evolutionary and molecular mecha- 
nisms underlying cross-species transmission is of critical 
importance for public health. After entering a new host 
cell, the success of a virus depends on its ability to 
overcome antiviral factors in the cell, such as protein 
kinase R (PKR). To investigate the process of virus 
transmission between species, we employed a recombi- 
nant vaccinia virus (WAEAK+RhTRSI) expressing the 
rhesus cytomegalovirus PKR antagonist RhTRSI. This 
protein inhibits some African green monkey (AGM) PKRs; 
however, it does not inhibit human or rhesus variants of 
PKR. Serial passaging WAEAK+RhTRSI in RhTRSI -resistant 
AGM cells resulted in rhtrsl duplication in the viral 
genome, which improved WAEAK+RhTRSI replication in 
AGM cells. Remarkably, rhtrsl duplication also enhanced 
virus replication in human and rhesus cells. In contrast, 
passage of WAEAK+RhTRSI in human cells, without prior 
adaptation in AGM cells, did not improve WAEAK+RhTRSI 
replication. These results support the hypothesis that 
amplification of a weak viral antagonist of a host defense 
protein in one species may enable cross-species transmis- 
sion into new hosts that are nonpermissive to the initial 
virus. 

by PKR, most virus families have evolved at least one mechanism 
to inhibit the PKR pathway [10]. 

This conflict between host restriction factors and their viral 
antagonists results in an "arms race" leading to rapid evolution of 
both sets of genes [11]. The extraordinarily high rate of positive 
selection (dN/dS) among primate PKR genes reflects the 
evolutionary pressure on PKR to evade virus antagonists 
[12,13]. To productively infect new host species, virus antagonists 
must rapidly adapt to escape these differences in PKR. In the 
current study, we modeled the process of viral adaptation through 



an intermediate, minimally permissive host. We experimentally 
evolved a recombinant vaccinia virus expressing the rhesus 
cytomegalovirus PKR antagonist rhtrsl in African green monkey 
cells expressing RhTRSI -resistant PKR. We demonstrate that 
amplification of the exogenous rhtrsl locus was an early adaption 
that is sufficient to rescue virus replication in minimally-permissive 
AGM cells. This amplification of the rhtrsl locus was also sufficient 
to expand the species tropism of these viruses, enabling them to 
infect both human and rhesus cells substantially better than the 
initial virus. Importantly, rhtrsl amplification did not occur when 
the initial virus was passaged through human fibroblasts, 
suggesting that amplification in AGM cells was a necessary 
intermediate step to expand the virus host range. Gene 
amplification is a universal mechanism of rapid adaptation 
occurring in eukaryotes [14,15], prokaryotes [16], and viruses 
[17,18], enabling diverse adaptations including evasion of host 
restriction factors (Reviewed in [19]). Our results suggest that gene 
amplification in an intermediate host may be a risk factor for 
broad cross-species transmission independent of other adaptive 
events. 

Results 

Variable inhibition of African green monkey PKR by 
RhCMV TRS1 

The rhesus cytomegalovirus (RhCMV) PKR inhibitor TRS1 
(RhTRSI) can block PKR activation and rescue replication of a 
vaccinia virus mutant lacking the PKR inhibitor E3L 
(WAE3L+RhTRSl) in several African green monkey (Chlorocehus 
aethiops, AGM) cell lines [20]. However, we discovered that a 
recombinant vaccinia virus expressing RhTRS 1 and lacking both 
of the known vaccinia PKR antagonists, E3L and K3L, 
(WAEAK+RhTRSI) produced 100 to 1000-fold less virus in 
AGM-derived PRO 11 90 cells relative to W-fJg (which contains 
both E3L and K3L) although it replicated almost as efficiently as 
W-pg in AGM-derived BSC40 cells (Fig. 1A). Thus, RhTRSI 
varies in its ability to support WAEAK replication in different 
AGM cells. 
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Figure 1. PKR inhibits WAEAK+RhTRSI replication in PRO1190 cells. (A) Titers of the indicated viruses produced 48 hpi (MOI = 0.1) of 
PR01 190 (gray bars) or BSC40 cells (black bars), measured using HFF+TRS1. (B) Immunoblot of lysates of PR01 190 stably transduced with lentiviral 
vectors expressing control (left) or PKR-specific (right) shRNAs, probed with anti-PKR or anti-actin antibodies. (C) Titers of the indicated viruses 48 hpi 
(MOI = 0.1) in either PRO1190+ctrl shRNA (grey bars) or PRO1190+PKR shRNA (black bars), measured using BSC40 cells. Results are represented as 
means +1 STD. 

doi:1 0.1 371/journal.ppat.1 004002.g001 
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This WAEAK+RhTRS 1 replication defect in PRO 1190 cells 
may be due to incomplete inhibition of PKR in these cells by 
RhTRSl. To test this hypothesis, we generated PRO 1190 cells 
stably expressing either a PKR-specific shRNA (PRO1190-PKR 
kd), which resulted in a 56% reduction of PKR expression, or a 
scrambled control shRNA (PROl 190-ctrl kd) (Fig. IB). Similarly, 
PKR-specific RT-qPCR demonstrated a 60% reduction in PKR 
mRNA from PROl 190-PKR kd cells relative to PROl 190-ctrl kd 
cells (2578 copies or 6503 copies PKR/ng total RNA respectively) 
but little difference between PROl 190 and PROl 190-ctrl kd cell 
PKR levels (Fig. SI). In PROl 190-PKR kd cells, 
WAEAK+RhTRS 1 replication was almost completely rescued 
to W-fSg levels (Fig. 1C). We detected a similar increase in 
WAEAK+RhTRS 1 replication after transiently transfecting 
PROl 190 cells with a siRNA specific for PKR, but not a control 
siRNA (data not shown). Sequence analysis of PROl 190 PKR 
identified three non-synonymous mutations relative to a previously 
reported AGM PKR (GenBank # EU733254) that is sensitive to 
RhTRSl (Fig. S2, [20]). Interestingly, one of these single 
nucleotide variants is heterozygous in PRO 1 1 90 cells, and changes 
a residue (T577M) that is evolving under positive selection in 
primates [12]. Although additional studies will be needed to 
determine whether one or more of these PRO 1190 PKR 
polymorphisms is responsible for increased resistance to RhTRS 1 , 
the results shown in Fig. 1 demonstrate that the block to 
WAEAK+RhTRS 1 replication in PRO 11 90 cells is mediated 
by PKR. 

Rhtrsl amplification expands virus tropism to human and 
rhesus macaque cells 

To determine whether WAEAK+RhTRS 1 could adapt to 
overcome the PKR-mediated block to replication in PRO 1 1 90 cells, 
we utilized a system of experimental evolution. We infected 
PRO 11 90 cells with WAEAK+RhTRS 1 at a low multiplicity of 
infection (MOI = 0.1). 48 hours post infection (hpi) we lysed the 
infected cells, titered the resulting virus, and infected new PRO 1 1 90 
cells, repeating this cycle multiple times. After four passages we 
observed a 10- to 100-fold increase in viral replication that 
remained stable for at least three subsequent passages in each of 
three independent lineages (Fig. 2). We next performed a 
competition assay to assess the relative fitness of the passaged virus 
in comparison to the initial WAEAK+RhTRS 1 virus. We co- 
infected PROl 190 cells with either WAEAK+RhTRS 1 or W-A, 




Figure 2. Experimental evolution of WAEAK+RhTRS 1 increases 
replication fitness in PROl 190 cells. PROl 190 cells were initially 
infected with WAEAK+RhTRS1 (MOI=0.1). Virus was harvested at 
48 hpi, titered on BSC40 cells, and used to infect fresh PRO1190 cells 
and the process was repeated. Three independent infections resulted in 
an ~10-fold to 100-fold gain of replication fitness in PR01 190 that was 
evident by passage four. 
doi:10.1371/journal.ppat.1004002.g002 



both of which express eGFP, and the same competitor, WAE3L 
+RhTRSl, which expresses (3-gal (MOI = 0.1 for each virus). Virus 
produced 48 hpi was titered on permissive BSC40 cells and the 
specific progeny viruses were enumerated by detecting fi-gal 
(W AE3L+RhTRS 1 ) and eGFP (WAEAK+RhTRS 1 and W-A) 
expression in the plaques (Fig. S3). WAEAK+RhTRS 1 replicated 
~ 3-fold better than WAE3L+RhTRS 1, whereas W-A replicated 
290-fold better than WAE3L+RhTRS 1 , confirming that serial 
passage through PRO 1190 cells increased the fitness of passaged 
viruses approximately 100-fold relative to the initial WAEAK 
+RI1TRSI. 

Since the passaged viruses replicated more efficiently in a 
minimally permissive cell line, we investigated the ability of the 
passaged viruses to replicate in primary cells from more divergent 
primates. We have shown that RhTRSl does not inhibit human 
or rhesus PKR in the context of WAE3L [20] . To determine 
whether the adaptations that evolved during serial passage in 
PRO 1 1 90 cells affected the virus species tropism, we infected 
primary human foreskin fibroblasts (HFF) or primary rhesus 
fibroblasts (RF) with WAEAK+RhTRS 1, each of the three 
passaged virus pools or W-fSg at an MOI of 0.1 (Fig. 3). As 
expected, WAEAK+RhTRS 1 replicated poorly and W-fjg 
replicated efficiently in each cell type. Remarkably, all three 
passaged pools replicated between 1000- to 10,000-fold better than 
WAEAK+RhTRS 1 in both HFF (Fig. 3, center) and RF (Fig. 3, 
right). For each pfu of W-A, W-B and W-C used to infect the 
cells, 5.2, 1.9, and 6.3 pfu of progeny emerged from HFF and 
40.8, 7.3, and 3.7 emerged from RF, respectively, suggesting that 
these viruses were sufficiently well adapted to enable continuous 
propagation in these cells (see below). However, these virus pools 
still replicated 10- to 100-fold better in PROl 190 than in either 
human or rhesus cells (Fig. 3, left). Thus, adaptation of 
WAEAK+RhTRS 1 in minimally permissive AGM fibroblasts 
also provides a substantial replication benefit in human and rhesus 
cells expressing distantly related PKR proteins. To elucidate the 
underlying mechanism for this gain of fitness, we harvested DNA 
from passage 7 viruses for genetic analyses and then passaged the 
viruses once more in PRO 11 90 to generate viral stocks for 
biochemical and infectivity analyses. 

Gene amplification as a mechanism of rapid adaptation in 
vaccinia virus has been well documented [1 7, 18,2 1] . To determine 
whether gene amplification could account for the broadly 
improved replication of passaged WAEAK+RhTRS 1 we per- 
formed paired-end Illumina based deep sequencing (Short Read 
Archive #SRP033208). Based on read depth, we detected 
duplication of the rhtrsl locus in all three passage 7 virus pools 
but not in WAEAK+RhTRS 1 (Fig. 4A). Each of the passaged 
pools contained between 1.4 and 1.9 copies of rhtrsl per genome, 
although these numbers reflect averages of a heterogeneous 
population of viral genomes. Confirming this estimate of rhtrsl 
copy number, the frequency of reads in which we detected a 
recombination site near the rhtrsl locus increased as a percentage 
of total reads in viruses predicted to have more copies of rhtrsl 
(Fig. 4B). 

We used PCR to confirm that the rhtrsl locus was amplified, 
using externally directed primers specific to rhtrsl that only amplify 
a product if there is a tandem duplication of the gene (Fig. S4A). 
We detected 3 kb products in all three virus pools, and 2.3 kb and 
1 .8 kb products only in the W-A virus pool (Fig. S4B). We were 
unable to obtain enough of the 2.3 kb band for further analysis, 
but we did characterize the 3 kb and 1 .8 kb products by Sanger 
sequencing. The larger product was identical in all three passaged 
virus pools, and represents a recombination between the vaccinia 
virus gene L5R upstream of rhtrsl with J2R downstream of rhtrsl. 
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Figure 3. Adaptation in minimally-permissive PRO1190 cells expands the species tropism of VVAEAK+RhTRSI . PRO1190 (left), HFF 
(center), or RF (right) were infected with the indicated viruses (MOI=0.1). Lysates were harvested 48 hpi and titered on BSC40 cells. Data are 
represented as means +1 STD. 
doi:1 0.1 371 /journal.ppat.1 004002.g003 
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Figure 4. Rhtrsl copy number variation relative to VVAEAK+ 
RhTRSI. (A) Relative read depth of lllumina sequencing reads in the 
passaged virus pools normalized to VVAEAK+RhTRSI. The graph is 
centered on the rhtrsl locus (grey bar). (B) Copy number of RhTRSI 
(rescaled so that parental strain = 1) plotted against breakpoint read 
frequency (per million mapped reads) demonstrates agreement 
between each independent measure of RhTRSI duplication. 
doi:1 0.1 371 /journal.ppat.1 004002.g004 



In the smaller product J2R recombined with the neoR gene, which 
was introduced as a selection marker during construction of 
WAEAK+RhTRS 1 (Fig. S3C). These two sites represented the 
predominant recombination sites (85.5% and 2.8% respectively) 
identified by lllumina deep sequencing. However, we found 
additional minor recombination sites by lllumina deep sequencing, 
including a 15 kb duplication in W-B. The presence of an 
identical recombination site in all three passaged virus pools 
suggests that duplication may have been present at a very low 
frequency in the initial virus population even though we did not 
detect it in the lllumina sequencing data. Regardless, taken 
together these data demonstrate that the copy number of rhtrsl in 
the viral genome was substantially increased by passage through 
PROl 190 cells. 

Unlike a previous study which identified adaptive point 
mutations arising after locus expansion [17], we did not detect 
any point or indel mutations in rhtrsl in any of the passaged virus 
pools. However, we identified 13 vaccinia virus gene mutations 
present at greater than 5 % frequency in at least one of the pools 
(Table 1). All three passaged pools had one or more of four 
different single nucleotide deletions within the A35R gene at 
frequencies ranging from 12 to 42%. Transition mutations 
affecting the A24R and A37R genes were present at >50% 
frequencies in W-A and W-B respectively, but were rare or 
absent in the other viral pools. None of the other mutations were 
detected in all three pools or occurred at >50% frequency in any 
pool, so are unlikely to account for the improved replication of the 
passaged viruses. However, the presence of these W gene 
mutations raised the question of whether the expanded species 
tropism we observed was due to the W gene mutations or to rhtrsl 
amplification. 

If rhtrsl amplification alone is sufficient for the observed increase 
in fitness, we reasoned that overexpression of rhtrsl in trans might 
rescue WAEAK+RhTRS 1 replication. To investigate this possi- 
bility, we stably transduced rhtrsl into HFF (HFF+RhTRSl), and 
confirmed RhTRSI expression by immunoblot (data not shown). 
We also prepared a control cell line (HFF-LHCX) by transducing 
the empty vector, LHCX, into HFF. In the control cells, 
WAEAK+RhTRS 1 replicated approximately 1000-fold less effi- 
ciently than W-Pg (Fig. 5A). In HFF+RhTRSl, WAEAK+ 
RhTRSI replication increased more than 100-fold. Thus, com- 
bined expression of RhTRSI from genes in both the cell and 
the infecting virus potentiated WAEAK+RhTRS 1 replication, 
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supporting the hypothesis that rhtrsl amplification alone is sufficient 
to expand the species tropism of WAEAK+RhTRS 1 . 

Although rhtrsl amplification provided a substantial growth 
benefit in HFF and RF, the passaged viruses still replicated 100- to 
1000-fold less efficiently than W-(5g (Fig. 3). To determine 
whether this incomplete rescue in HFF was due to incomplete 
PKR inhibition or represented a second block to replication, we 
infected HFF stably transduced with either a PKR specific shRNA 
(HFF-PKR kd) that reduces PKR expression >95%, or a non- 
specific shRNA (HFF-ctrl kd) [22]. Knocking down PKR 
increased WAEAK+RhTRS 1 replication -1000-fold, indicating 
that PKR is a major barrier to replication in these cells (Fig. 5B). 
All three passaged virus pools replicated ~ 10-fold better in HFF- 
PKR kd cells than in HFF-ctrl kd cells, suggesting that rhtrsl 
amplification, which fully inhibits PRO 1190 PKR, only partially 
inhibits human PKR. However, these viruses all replicated ~ 1 0- 
fold less well than W-pg in the HFF-PKR kd cells. This 
remaining replication defect may be due to incomplete PKR 
knockdown in these cells [20], although it is also possible that an 
additional host factor inhibits WAEAK+RhTRS 1 replication in 
HFF. 

RhTRSI inhibits the PKR pathway after 
autophosphorylation in multiple species 

Unlike many known PKR inhibitors that block PKR phos- 
phorylation [10,23], RhTRSI inhibits the PKR pathway after 
PKR phosphorylation but prior to eIF2a phosphorylation [20]. 
However, it is possible that RhTRSI amplification inhibits PKR 
through an alternative mechanism, such as dsRNA sequestration. 
To determine whether rhtrsl amplification altered the mechanism 
of PKR inhibition, we infected PRO 1 190 with W-pg, WAE.3L, 
WAEAK+RhTRS 1, and W-A at an MOI of 3, and coUected 
cell lysates 24 hpi. 35 S metabolic labeling demonstrated that 
WAEAK+RhTRS 1 expressed vaccinia virus proteins in 
PRO 1190 cells, though in much lower quantities than W-(3g, 
consistent with the former virus being unable to inhibit PKR 
completely (Fig. 6, top left panel). In contrast, W-A produced 
abundant vaccinia virus proteins, similar to W-(3g, confirming 
that this virus efficiendy inhibits the PKR pathway. 

We used immunoblot analyses to determine the stage of the 
PKR pathway inhibited by each virus (Fig. 6, lower panels). 
RhTRSI expression was noticeably higher in the W-A infected 
cells compared to those infected with WAEAK+RhTRS 1 , 
consistent with the better replication of W-A in PRO 1 1 90 cells 
(Fig. 6, lanes 4 and 5). PKR phosphorylation was elevated in all 
virus infected cells except W-|3g, consistent with our previous 
report [20] that even when RhTRSI blocks the PKR pathway, it 
does not block PKR autophosphorylation (Fig. 6, lanes 3-5). 
Phospho-eIF2a levels in WAEAK+RhTRS 1 infected cells were 
intermediate between mock and WAE3L infected PRO 1190 
cells, suggesting that a single RhTRSI gene weakly inhibits the 
PKR pathway (Fig. 6, lane 4). Infection with W-A resulted in low 
levels of eIF2ot phosphorylation, similar to that detected in W-fig 
infected PROl 190 cells (Fig. 6, compare lanes 2 and 5), indicating 
that RhTRSI amplification is sufficient to completely inhibit 
PKR-mediated translational shutdown at a stage after PKR 
phosphorylation. Together, these data suggest that PKR-mediated 
pressure in PRO 1190 cells selected for rapid amplification of the 
rhtrsl locus, and that this amplification was sufficient to enable the 
virus to block PKR-mediated defenses in PRO 1 1 90 cells without 
altering the mechanism of RhTRS 1 -mediated PKR inhibition. 

In previous studies, RhTRSI alone was insufficient to inhibit 
human PKR [20]; however, the passaged viruses reported here 
replicate substantially better in HFF than WAEAK+RhTRS 1 , 
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Figure 5. RhTRSI overexpression or PKR knockdown permits efficient VVAEAK+RHTRS1 replication in HFF. (A) Viral production 48 hpi 
of control (HFF+LHCX, grey bars) or HFF expressing RhTRSI (HFF+RhTRS1, black bars) with the indicated viruses (MOI = 0.1), determined by titering on 
HFF+TRS1. (B) Viral production 48 hpi of control (HFF-ctrl kd) or HFF with PKR knocked down by shRNA (HFF-PKR kd) with the indicated viruses 
(MOI = 0.1), determined by titering on BSC40 cells. Data are represented as mean +1 STD. 
doi:1 0.1 371 /journal.ppat.1 004002.g005 



suggesting that amplification of rhtrsl is able to inhibit PKR at least 
partially in these cells. To test this hypothesis and determine 
whether the mechanism of PKR inhibition was the same in HFF 
as it is in PRO 1 190, we infected HFF (MOI = 3) and prepared cell 
lysates 24 hpi. In contrast to infection of PRO 1 190 cells, infection 
of HFF with WAEAK+RhTRS 1 resulted in nearly complete shut 
off of protein synthesis by 24 hpi (Fig. 6, top right panel) and 
produced only trace amounts of RhTRSI. Compared to 
WAEAK+RhTRS 1, W-A infection of HFF resulted in detect- 
able, though still low overall levels of ' 5 S labeled proteins, and 
much more RhTRSI (Fig. 6, lanes 9 and 10). PKR phosphory- 
lation and eIF2ct phosphorylation were elevated after infection 
with both WAEAK+RhTRS 1 and W-A compared to mock or 
W-Pg controls (Fig. 6, compare lane 9 to lanes 6 and 7). These 
data suggest that a single copy of rhtrsl was insufficient to inhibit 
the translational shutoff mediated by human PKR, but amplifi- 
cation of this weak antagonist resulted in partial inhibition of 
human PKR allowing enough protein synthesis to support a 
modest level of virus replication. 

RhTRSI amplification in PRO1190 is a necessary 
intermediate for WAEAK+RhTRS1 replication in HFF 

Because selection in AGM cells resulted in a broad expansion of 
viral species tropism, we investigated whether passage of 
WAEAK+RhTRS 1 direcdy in HFF would similarly select for 
mutants, such as rhtrsl amplification, that improved replication in 
HFF. We therefore seriaUy infected PROl 190 cells and HFF with 
WAEAK+RhTRS 1 in parallel (Fig. 7). In PRO 11 90 cells, viral 
fitness again increased after four passages. In contrast, virus 
replication was strongly inhibited in HFF, and we were unable to 
detect any viral replication after three rounds of infection in all 
three pools. These data suggest that, under these experimental 
conditions, adaptation in PROl 190 cells was a necessary 
intermediate step for improved replication in HFF. 

Finally, we evaluated whether the improved replication of the 
viruses that had been passaged in PRO 1 1 90 cells was sufficient to 
enable stable propagation in HFF. We therefore serially infected 
HFF with W-A, W-B and W-C at low multiplicity of infection 
(MOI = 0.1 at each passage). We were able to propagate these 
viruses in HFF for at least four passages. Moreover, replication 
increased between 5- to 14-fold after only two passages, suggesting 
further adaptation occurred in HFF cells (Fig. 7B). To define 
mutations that may have evolved during serial passage in HFF, we 
performed paired-end Illumina based deep sequencing on viral 



DNA isolated after the fourth round of passage in HFF. Again, we 
did not find any mutations in rhtrsl; however, we detected an 
average rhtrsl copy number of 1.7, 2.8, and 2.7 for the virus pools 
derived from W-A, W-B and W-C, respectively, representing 
an average expansion of the rhtrsl locus by approximately one 
additional copy relative to the PRO 1 1 90-adapted viruses. In 
addition, two W gene mutations (in F7L [indel] and J6R) that 
arose during PRO 1190 adaptation were lost after HFF adaptation, 
and two new mutations (in F7L [missense] and H4L) of uncertain 
significance appeared during HFF adaptation (Table S3). Taken 
together, our study suggests the hypothesis that gene amplification 
acts broadly to increase replication in a variety of hosts, and may 
provide a molecular foothold that allows for continued species- 
specific adaptation of the virus in more resistant hosts. 

Discussion 

Cross-species pathogen transmissions have been responsible for 
more than 60% of all emerging infectious diseases in humans 
during the past 70 years [24]. Human immunodeficiency viruses, 
avian influenza viruses, and the recently described Middle East 
respiratory syndrome coronavirus exemplify the ongoing threat 
and potential of animal viruses to spread to and among humans, 
highlighting the urgent need to understand the mechanisms 
underlying cross-species transmission and adaptation to new hosts. 
One such mechanism, genetic locus amplification in response to 
selective pressure, has been observed in both viruses and bacteria 
[16-18,21]. Here we demonstrate that amplification of the 
exogenous gene rhtrsl is sufficient to block potent PKR-mediated 
inhibition and improve WAEAK+RhTRS 1 replication in AGM- 
derived PROl 190 cells. This adaptation also expanded the species 
tropism of the virus, enabling markedly improved replication in 
otherwise resistant human and rhesus monkey cells. Importantly, 
WAEAK+RhTRS 1 failed to replicate in HFF to a level sufficient 
to sustain transmission upon serial passage, demonstrating that 
adaptation in PRO 1 1 90 was a critical intermediate step to expand 
the viral species tropism. Thus, the process of adaptation in one 
host may increase the likelihood of virus transmission to a variety 
of divergent species. 

Under PKR-mediated selective pressure, the rhtrsl locus 
amplified during serial passage of WAEAK+RhTRS 1 in 
minimally permissive PRO 11 90 cells (Fig. 4). It is not clear 
whether the initial duplication(s) occurred during preparation of 
the WAEAK+RhTRS 1 stock in BSC40 cells or during the first 
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Figure 6. Effects of RHTRS1 overexpression on the PKR pathway in PRO1 190 or HFF. PRO1 190 (left) or HFF (right) were infected with the 
indicated viruses (MOI = 3) for 24 hpi. Proteins were labeled with 35 S-methionine for 1 h and cell lysates were analyzed by autoradiography (top 
panel) or immunoblotting with the indicated antibodies (bottom panels). 
doi:1 0.1 371 /journal.ppat.1 004002.g006 



few passages in the PRO 1190 cells. The observation of a faint 
PCR product from the starting virus, VVAEAK+RhTRS 1 using 
outward directed rhtrsl primers in one of three experiments (Fig. 
S4), and the detection of an identical recombination break point in 
all 3 independently passaged virus pools suggest the amplification 
may have been present at a very low level in the starting virus. 
However, amplification of the locus in WAEAK+RhTRSl 
occurred below the level of Illumina deep sequencing detection 
(Fig. 4), supporting the idea that if rhtrsl duplications are present 
they are rare in the initial WAEAK+RhTRSl stock. Addition- 
ally, we detected recombination between J2R and neoR (Fig. S4) 
only in the W-A pool, suggesting that recombination events did 
arise during serial passage. Regardless of when they arose, 
amplifications of the rhtrsl locus were substantially enriched 
during virus passage under selective pressure (Figs. 2 and 4). A 
previous study demonstrated amplification of K3L as an adaptive 
mechanism against human PKR with strikingly similar kinetics to 
the current study [17]. Together, these two studies support the 
hypothesis that preexisting or frequently arising de novo gene 
duplications enable vaccinia virus to adapt rapidly to selective 
conditions imposed by relatively resistant host restriction factors. 

The "accordion hypothesis" of rapid evolution posits that gene 
amplification provides a replication benefit to the virus and those 
extra copies of a weak viral antagonist of host defenses provide 
additional templates to acquire potentially adaptive mutations. 



Indeed, Elde, et al. detected such an adaptive mutation in K3L 
(H47R), apparently arising after amplification of the locus [17]. 
Therefore, we were surprised that no mutations arose in rhtrsl, 
although it may be that additional rounds of replication would 
reveal such mutations. We did, however, detect 13 mutations in 
endogenous vaccinia virus genes that occurred at >5% frequency 
(Fig. 4). While our data suggest that none of these mutations are 
necessary to expand the species tropism of WAEAK+RhTRS 1 
(Fig. 5), we have not ruled out the possibility that they may provide 
some replication benefit. The presence of mutations in A24R, 
A35R, and A37R are the most intriguing, as they were either 
present at >50% frequency in one pool (A24R and A37R) or 
detectable in all three pools (A35R). None of these genes has been 
previously implicated as a PKR antagonist. A37R is conserved 
across multiple poxvirus families (http://www.poxvirus.org), but 
its function is unknown. A24R is a subunit of RNA polymerase. If 
the mutation we identified acts like some other reported RNA 
polymerase mutations to decrease transcription elongation 
[25,26], this mutation might result in less dsRNA production 
and therefore less PKR activation. A35R is a gene of unknown 
biochemical function that may be involved in evasion of the 
adaptive immune response [27,28]. All three passaged virus pools 
contained nucleotide deletions in A35R at greater than 10% 
frequency. A35R orthologs are conserved across most poxviruses 
but, intriguingly, variola virus contains a truncation in its A35R 
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Figure 7. PROI 190-adapted viruses, but not VVAEAK+RHTRS1, stably replicate in HFF. (A) PRO1190 or HFF were infected with 
VVAEAK+RhTRSI (MOI = 0.1), for 48 h, after which viral progeny were collected, titered on BSC40 cells, and used to infect new cells (MOI = 0.1, or the 
full volume of the harvested lysate when too little virus was obtained). Three independent infections resulted in a >1 0-fold gain of replication fitness 
in PRO1 190 over the course of the experiment, but no virus survived three rounds of passage through HFF. (B) HFF were infected with PR01 190- 
adapted VV-A, W-B, or VV-C virus pools (MOI = 0.1). Virus progeny harvested at 48 hpi was titered on BSC40 cells and used to infect new HFF and the 
process was repeated. All three virus pools replicated in HFF for at least four rounds of serial passage, and acquired an additional, modest gain in 
replication fitness. 

doi:1 0.1 371 /journal.ppat.1 004002.g007 



gene, demonstrating that similar truncations have evolved in the 
past. Further studies are underway to evaluate the potential 
contributions of these mutations to viral replication. 

Consistent with previous studies of RhTRSl in other AGM cell 
types [20], PKR, but not eIF2ot, was phosphorylated in PROI 190 
cells infected with W-A. Thus, the improved replication of the 
passaged virus pools in PRO 1 1 90 cells was likely due to 
enhancement of the basic RhTRS 1 -mediated inhibition of PKR, 
and not due to another mechanism, such as dsRNA sequestration 
or reduction in the abundance of dsRNA. WAEAK+RhTRS 1 
infected PROI 190 cells had phospho-eIF2ot levels lower than that 
found in WAE3L infected cells (Fig. 6, lanes 3 and 4), suggesting 
that even just a single copy of rhtrsl is able to inhibit PKR function 
to a small degree, but amplification of the locus appears to be 
needed to express enough RhTRSl to inhibit eIF2ot phosphor- 
ylation potendy and enable efficient viral replication. It is also 
possible that elevated expression of RhTRS 1 is necessary to block 
another activity of PKR, such as autophagy or inflammasome 
responses [29-31] that might aid in replication. 

Adaptation of WAEAK+RhTRS 1 to PROI 190 cells provided 
a substantial replication advantage in both human and rhesus 
monkey fibroblasts. Although W-A replicated to much higher 
titers and expressed more RhTRSl than WAEAK+RhTRS 1 in 
HFF, 35 S metabolic labeling revealed relatively low protein 
synthesis rates and eIF2ot phosphorylation was still elevated after 
W-A infection compared to W-fig. A substantial proportion of 
the block to W-A replication in HFF is still mediated by PKR 
despite RhTRSl overexpression (Fig. 5B), suggesting that further 
adaptation in HFF may be necessary to block the PKR pathway in 
HFF completely. These results indicate that RhTRS 1 overexpres- 
sion blocks human PKR incompletely. Nonetheless, viruses that 
had adapted by initial passage in PROI 190s replicated in HFF at 
a level sufficient to enable sustained passage in HFF (Fig.7B). 
Furthermore, these PRO 1 1 90-adapted viruses acquired additional 
changes as a result of serial passage in HFF, although the 
biological relevance of these changes is currently unclear. 



Combined with the observation that serial passage of 
WAEAK+RhTRS 1 in HFF failed to generate any adapted 
viruses (Fig. 7A), these data suggest that adaptation of rhtrsl in 
minimally permissive PRO 1 1 90 cells was a critical intermediate 
step in the generation of a virus with broadened species tropism. 

Cross-species pathogen transmission is an important source of 
emerging infections worldwide. Our study illustrates that gene 
amplification of a weak viral antagonist of PKR can broaden the 
host range of vaccinia virus. The presence of gene families in other 
large DNA viruses, e.g. the cytomegalovirus US22 family, of which 
rhtrsl is a member, provides indirect evidence that episodes of 
locus amplification have also occurred in other viruses. Adaptation 
of a viral antagonist through non-synonymous mutations has the 
potential to confer a species-specific advantage for the virus in a 
specific host. In contrast, gene amplification and subsequent over- 
expression of the antagonist is more likely to increase protein 
activity through mass action effects in a variety of hosts. Thus, 
gene amplification may be a common evolutionary strategy 
employed by large DNA viruses, permitting modest replication 
in otherwise resistant host species and providing a molecular 
foothold that enables further adaptations necessary for more 
efficient replication and spread in the new host. 

Materials and Methods 

Viruses and cell culture 

AGM fibroblasts (PRO 11 90, Coriell Institute for Medical 
Research) human foreskin fibroblasts (HFF), rhesus fibroblasts 
(RF), BSC40 and HeLa cells, and derivative cell lines were 
maintained in Dulbecco's modified Eagle's medium supplemented 
with 10% NuSerum (BD Biosciences) as previously described [32]. 
At times, PRO 1 1 90 cell lines were also propagated in Minimal 
Essential Medium with 20% fetal calf serum and anitbiotics to 
enable more rapid growth. PRO 1 1 90 cells were transduced with a 
nonsilencing control or PKR-targeting shRNA lentiviral vectors 
(Open Biosystems, catalogue numbers RHS4430-988 19555 and 
RHS4346, respectively) and selected in puromycin (5 ug/ml) to 
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generate PRO1190-ctrl k/d and PRO1190-PKR k/d cells. 
pEOJ364 was constructed by moving the RhTRSl gene, with a 
C-terminal biotinylation signal and 6x-His tag, as a Hindlll/Pmel 
fragment from pEQ1215 [20] into the Hindlll/Hpal sites of 
pLHCX (Clontech Laboratories, Inc). HFF-LHCX and 
HFF+RhTRS 1 were produced by transducing HFF with retroviral 
vectors made using LHCX and pEQJ364, respectively and 
selecting with hygromycin B (100 ug/ml). 

Vaccinia virus (W) Copenhagen strain (VC2) [33] and 
WAE3L [34], both obtained from Bertram Jacobs (Arizona 
State University), and W-fSg (VC2-LacZ in [20] were propagated 
and titered in BSC40 cells. 

VC-R2 (W AE3LAK3L) was constructed by replacing the E3L 
gene in the K3L-deleted VACV vP872 strain (AK3L in VC2 
background, provided by Bertram Jacobs) [35] by homologous 
recombination. The 518 bp 5' arm was created by PCR 
amplification of VC2 DNA with primers C15 (5 ' -GATTAAGGG- 
TACTAGC GGC ACC G ' 3 ') x C 1 6 (5' -TTTT AG AG AG AAC TA- 
ACACAACCAGC-3'). The 512 bp 3' arm was created by PCR 
amplification of VC2 DNA with primers C19 (5'-GTGTAG- 
TAAGCTAGCGAGCTCGGTACCTTCTAGTTATCAATAA- 
C AGTT AGT AGTTT AG- 3 ' ) xC20 (5'-CCAACAAACTGTTC- 
TCTTATGAATCG-3'). The reading frame of EGFP including 
the PEST sequence was amplified with primers C17 (5'-GCT- 
GGTTGTGTTAGTTCTCTCTAAAACCCGGGATCCACCG- 
GTCGCC-3')xC18 (5'-GGTACCGAGCTCGCTAGCTTAC- 
T AC AC ATTG ATC C T AGC AG AAGC - 3 ' ) using pD2EGFP-Nl 
(Clonetech) as the template. PCR products were gel-purified and 
mixed together as template for fusion PCR using C15xC20. 
PfuUltra polymerase (Agilent Technologies) was used for these PCR 
reactions. PCR products were cloned into the pCR2.1 TOPO 
vector to generate plasmid S96. S96 was used as template for PCR 
amplification of marker +5' and 3' arms using C15xC20 followed 
by gel-purification. BS-C-1 cells grown on 12 well plates were 
infected with vP872 at MOI = 2 and transfected 2 hours after 
infection with 1 |j,g of the purified PCR product. Cell lysates were 
collected after 18 hours and plated in a dilutions series on 
RK13+E3L+K3L cells [36]. Green plaques were picked after 
48 hours at the highest dilution possible and plaque purified an 
additional three times on RK13+E3L+K3L cells. EGFP expression 
in VC-R2 is under the control of the endogenous E3L promoter. 
WAEAK was propagated and titered using HFF+TRS1 cells (HF- 
TRS1 in [37]). 

WAEAK+RhTRS 1 was constructed by homologous recombi- 
nation of plasmid pEQJ233 [20] into the thymidine kinase (TK) 
locus of WAEAK. Recombinant virus was plaque purified three 
times in BSC40 under G418 selection and subsequently propa- 
gated and titered on BSC40 cells. 

PRO1190 PKR sequence analysis 

Total RNA from PRO 1190 cells was amplified using previously 
reported PKR specific primers [12]. The amplification product was 
gel purified and cloned using the StrataClone PCR cloning kit 
(Agilent). Multiple plasmids were submitted for Sanger sequencing 
using PKR specific sequencing primers (#859: 5'-ATGGCTGGT- 
GATCTTGCAC; #860: 5 ' -GTGAAC AACTC ACTTGCTTC ; 
#861: 5 ' -GAAACTAGAC AAAGTTTTGGC ; #862: 5'-CTAA- 
CATGTATGTCGTTCCT; #863: 5'-AAGGCACTTAGTCTT- 
TGATC; #864: 5'-TCTGATATCTCAAGCAATGC). Contigs 
were assembled and curated in Geneious Pro v4.8.5 (GenBank 
#KF728076-7). The predicted amino acid sequence was aligned to 
the predicted amino acid sequences of previously reported AGM 
(GenBank # EU733254), rhesus macaque (GenBank# EU733261), 



and human PKR (GenBank # NM001 135651) using ClustalW2 
(http:/ / www.ch.embnet.org/software/ClustalW.html) [38] . 

Experimental evolution of WAEAK+RhTRS1 

For every round of infection, triplicate confluent 1 0 cm dishes 
(Figure 2) or 6-well plates (Figure 7) of PROl 190 or HFF were 
infected with WAEAK+RhTRS 1 (MOI = 0.1). Two days post- 
infection, cells were collected, pelleted and resuspended in 1 mL 
DMEM+10% NuSerum. After three freeze/thaw cycles, virus 
titers were determined on BSC40 by plaque assays and used for 
the next round of infection. If insufficient virus was produced to 
infect cells at MOI = 0.1, the entire volume of virus lysate was used 
in the subsequent round of infection. 

Vaccinia virus DNA from passage 7 in PRO 1 1 90 (Figure 2) or 
passage 4 in HFF (Figure 7 right) viruses was purified from 
infected cell cytoplasmic extracts for genetic analyses described 
below [39]. The passage 7 pools were further expanded (passage 8) 
in PROl 190 for use in virologic assays. 

Competition assay 

PRO 11 90 cells were co-infected with 0.1 MOI of either 
WAEAK+RhTRS 1 or W-A, and 0.1 MOI of WAE3L+ 
RhTRSl as a common competitor. Two days post-infection cells 
were collected, pelleted and resuspended in 1 mL DMEM+10% 
NuSerum. After three freeze/thaw cycles, virus titers were 
determined on BSC40 by plaque assays. WAE3L+RhTRS 1 
plaques were detected with the fS-gal substrate ImaGene Red 
Ci 2 RG following the manufacturer's directions (Life Technologies). 
Plaques were imaged on a Typhoon Trio imager (eGFP - 488 nm 
excitation, 520 BP 40 filter; ImaGene Red - 532 nm excitation, 580 
BP 30 filter) at 50 |J,m/pixel resolution, and classified as GFP+, 
ImaGene Red+, or double positive using ImageJ software (http:// 
rsb.info.nih.gov/ij/). 

Genomic analyses 

Externally directed PCR. Externally directed oligonucleo- 
tides were designed that bind near the 5' or 3' end of RhTRSl 
(#915: 5 ' -TGTGGGAGGATGC ATTGC AG; #916: 5'-GGC- 
GACTACAATCCCCATTG). High-fidelity PCR was performed 
with 1 00 ng of purified virus DNA following the manufacturer's 
directions (Phusion, Thermo Scientific). Reactions were pulsed at 
98°Cx30s followed by 25 cycles of 98°Cxl0s, 65°Cx20s, 
72°C x3 min. Amplification products were gel-purified and cloned 
using the Strataclone PCR cloning kit (Agilent) following the 
manufacturer's suggested protocol. Recombination sites were 
identified by Sanger sequencing of these products. 

Library preparation and sequencing. From each viral 
pool, 75 ng of viral DNA was sheared and ligated to adaptors by 
transposition using the Nextera kit (Epicentre), following the 
manufacturer's direction with several modifications: transposition 
was performed in 20 ul volume using 0.1 ul Nextera transposase 
enzyme and "HMW" reaction buffer at final IX. Transposition 
was carried out for 10 minutes at 55°C, after which sheared, 
adaptor-ligated templates were transferred directly to PCR 
reactions for amplification. In addition, PCR incorporated a 
sample-specific barcode tag on the reverse primer. Amplified 
libraries were pooled and cleaned using Ampure beads (Beckman 
Coulter) at a 1.8:1 ratio of beads to input. The pooled libraries 
were sequenced on an Illumina HiSeq 2000 instrument (Illumina) 
using 101-bp forward and reverse reads with a 9-bp index to read 
the per-sample barcode. 

Sequence analysis. De novo assembly guided by the 
Vaccinia virus Copenhagen genome (GenBank #M35027.1) was 
used to construct a reference assembly sequence for 
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WAEAK+RhTRS 1 (Tables SI and S2). Shotgun reads from the 
parental virus were assembled to contigs using ABySS 1.3.4 [40]; 
parental virus reads were then realigned against the resulting 
assembly contigs using bwa 0.6.1 [41]. Discrepancies between 
reads and contigs, reflecting de novo assembly errors, were corrected 
in the contigs following manual examination, and contigs were 
ordered and concatenated by comparison to the Copenhagen 
sequence. 

Reads from each evolved viral pool were aligned to the parental 
reference genome, and single-base and short-indel variants were 
called using the Genome Analysis Toolkit [42]. Custom scripts 
were used to determine the allele frequency (fraction reads 
supporting variant allele out of all aligned reads at a given site) for 
all variants found in any sample. Copy number was estimated by 
counting the depth of read coverage within sliding windows of 
1 00 bp, correcting for effects of differences in G+C composition, 
and dividing by the parental strain read depth. Everted reads at 
duplication breakpoints were identified using bwasw [43] and 
custom scripts, as previously described [17]. 

Immunoblot analyses 

Cells were mock infected or infected with vaccinia viruses 
(MOI = 3). One day postinfection, the cells were lysed in 2% 
sodium dodecyl sulfate (SDS). Equivalent amounts of the lysates 
were separated on 10% SDS-polyacrylamide gels, transferred to 
polyvinylidene difluoride (PVDF) membranes, and probed with 
one of the following antibodies: PKR (sc-6282; Santa Cruz 
Biotechnology, Inc.), phospho-PKR (T446; 1120-1; Epitomics), 
eIF2ot or phospho-eIF2a (Ser51) antibody (both from Cell 
Signaling Technology, catalog numbers 9722 and 9721, respec- 
tively), TRS1 a999 [37], or actin (A2066; Sigma). All purchased 
antibodies were used according to the manufacturer's recommen- 
dations. Proteins were detected using the Western Star chemilu- 
minescent detection system (Applied Biosystems) according to the 
manufacturer's recommendations. Densitometry measurements 
were performed using ImageJ. 

RT-qPCR 

Total RNA was extracted from PRO 11 90, PRO1190-ctrl kd, 
and PRO1190-PKR kd cells with TRIzol reagent following the 
manufacturer's protocol (Invitrogen), and 1 ng of total RNA was 
assayed per reaction. The standard curve was generated from 10- 
fold serial dilutions of a PRO 1 1 90 PKR containing plasmid 
(described above in PRO1190 PKR sequence analysis, pEQ_1334) 
diluted in 200 pg/uL salmon sperm DNA. All samples were 
amplified in triplicate using the GoTaq 1-step RT-qPCR kit 
following the manufacturer's protocol (Promega) using PKR 
specific primers (#1063: 5 ' -C AC AGAATTGACGGAAAGAC ; 
#1064: 5'-ATCCCAACAGCCATTGTAGT). RT-qPCR was 
performed on a Rotor-Gene Q, thermocycler (Qiagen) with 
temperature holds at 37°Cxl5 min and 95°Cxl0 min followed 
by 40 cycles of 95°CxlO s, 60°Cx30 s. Raw data was analyzed 
using the included Rotor-Gene Q_ series software using the 
automatic cycle threshold (Ct) setting for assigning baseline and 
threshold for Ct determination. Nonlinear regression analysis to 
determine PKR copy number in the experimental samples was 
performed in GraphPad Prism 6. 

Metabolic labeling 

PRO 1 1 90 and HFF cells were mock infected or infected with 
the indicated viruses (MOI = 3). At 24 hpi, the cells were labeled 
for 1 h with 100 uCi/ml L-[ 35 S]methionine/L-[ 35 S]cysteine 
(EasyTag express protein labeling mix; PerkinElmer) in medium 
lacking methionine and cysteine. The cells were then lysed in 2% 



SDS. Equivalent amounts of protein from each sample were 
separated on 10% SDS-polyacrylamide gels, dried, and visualized 
by autoradiography. 

Supporting Information 

Figure SI RT-qPCR analysis of PKR copy number in 
PRO1190-derived cell lines. RT-qPCR was performed on 
serial 10-fold dilutions of a PKR containing plasmid from 10 to 
10 1 copies of PKR to generate a standard curve (top). Data are 
represented as the mean +/— the 95% confidence interval. PKR 
copy number/ng total RNA was interpolated from this standard 
curve for PRO 1190, PRO1190-ctrl kd, and PRO1190-PKR kd 
cells (bottom). Relative to PRO 11 90 cells, PRO1190-ctrl kd cells 
expressed slightly less PKR (20% decrease), while PROl 190-PKR 
kd cells had a 69% decrease in PKR expression. Similarly, 
PROl 190-PKR kd cells expressed 60% less PKR than PKR-ctrl 
kd cells, consistent with Fig. IB. 
(TIF) 

Figure S2 Amino acid alignment of African green 
monkey, human, and rhesus PKR. Predicted PKR amino 
acid sequences from AGM (PROl 190T1610 and EU733254), 
rhesus macaque (GenBank# EU733261), and human (GenBank 
# NM001 135651) were aligned using CLUSTAL W (1.83) [38]. 
We identified three non-synonymous amino acid differences (grey 
or green highlighted residues) between the two AGM alleles. One 
of these differences is at a site that is evolving under positive 
selection in primates (green highlighted residue) [12]. Relative to 
PROl 190 PKR, rhesus and human PKR are 90.8% and 83.2% 
identical at the amino acid level, respectively. 
(TIF) 

Figure S3 Increased fitness of VV-A compared to 
WAEAK+RhTRSl assessed by indirect competition 
assay. PRO 1 1 90 cells were co-infected with either 
WAEAK+RhTRSl or W-A (MOI = 0.1) and the same 
competitor virus, WAE3L+RhTRS 1 (MOI = 0.1). Two days 
post-infection viral progeny were collected and titered on 
BSC40 cells. (A) Representative Typhoon image of virus titer 
plates produced by infection with WAE+RhTRSl in combi- 
nation with either WAEAK+RhTRSl (top panels) or W- 
A (bottom panels). (B) Plaques were scored for eGFP expression 
(WAEAK+RhTRSl or W-A), [5-gal expression (WAE3L+ 
RhTRSl) or both (double +) as described in Materials and 
Methods. W-A replicated ~ 100-fold better than WAEAK+ 
RhTRSl relative to WAE3L+RhTRS 1 . Data are represented as 
the mean +1 SD. 
(TIF) 

Figure S4 Predominant recombination sites identified 
in rhtrsl locus amplification. (A) Schematic of externally 
directed PCR. Externally directed oligonucleotides (grey arrows) 
were designed to bind to the ends of rhtrsl (white). Only tandem 
duplications will produce an amplification product. (B) Externally 
directed PCR revealed enrichment of rhtrsl amplification products 
in the passaged virus pools. (C) Schematic of the initial rhtrsl locus 
in WAEAK+RhTRS 1 (top), and the predominant recombination 
sites identified in all three passaged viruses (middle), or in W-A 
alone (bottom). Vaccinia virus genes are colored green or blue, 
and exogenous sequences are colored dark red (rhtrsl) or brown 
(neoR). PCR products amplified by externally directed primers 
(black bars), the duplicated sequence (grey boxes), and the 
recombination sites (red arrowheads), are indicated. Recombina- 
tion sites occur (middle) after nt 84267 (J2R) and before nt 83200 
(L5R) relative to the reference vaccinia virus. Copenhagen strain 
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(Genbank #M35027.1), or (bottom) after nt 84237 (J2R) and 

before nt 622 (neoR - pcDNA3.1, Invitrogen). 

(TIF) 

Table SI Genomic coordinates for regions in whole 
genome alignment between Vaccinia Copenhagen refer- 
ence (GenBank M35027.1) and parental genome 
(WAEAK+RhTRSl genome) sequenced here, as deter- 
mined by MUMmer 3.2.1 [44]. 
(DOCX) 

Table S2 Substitution and short indel variants (<10 bp) 
present in parental genome (WAEAK+RhTRSl) as 
compared to Vaccinia Copenhagen genome (excluding 
any in the inverted terminal repeat regions). Predicted 
effects upon protein coding genes as annotated in the 
Copenhagen reference are listed, with more than one line 
per variant indicating effects upon multiple overlapping gene 
models. 
(DOCX) 
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